Positronium Spectroscopy 
in 

a Magnetic Field 



Jan Govaerts 

Institut de Physique Nucleaire 
Universite Catholique de Louvain 
2, Chemin du Cyclotron 
B-I348 Louvain-la-Neuve, Belgium 



Abstract 

Hyperfine spectroscopy of positronium formed in the presence of a static mag- 
netic field is considered. Generalising the situation hitherto developed in the 
literature, the magnetic field is not assumed to be parallel to the momentum 
of incoming polarised positrons, while the possibility of electron polarisation is 
also included in the analysis. The results are of application to high sensitivity 
positron polarimeters used in current (3 decay experiments. 
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1 Introduction 



Hyperfine structure in positronium is a basic physical fact with wide ranging imph- 
cations. The energy difference between singlet and triplet states, and their lifetimes 
dominated by disintegrations into two and three photons, respectively, are observ- 
ables providing crucial tests for quantum electrodynamics (QED)[|l| and the Standard 
Model [0 of electroweak and strong interactions. On a more practical level, positronium 
physics is also essential in a series of technological developments both in the design of 
detectors for particle physics experiments and in applications to problems of solid state 
physics (for a detailed review with references to the original literature, see Refs.|l|, |^). 

In particular, hyperfine structure of positronium has been put to use since over 
thirty years|^ in measurements of positron polarisations. One instance where such 
experiments come immediately to bear on the structure of fundamental interactions 
is in the case of positrons emitted in the [3 decay of radioactive nuclei. Specifically, 
any deviation from purely right-handed polarisation of positrons emitted in the j3 
decay of polarised nuclei would point to physics beyond the Standard Model in its 
electroweak sector. Actually, experiments with this purpose in mind and using such 
a positronium based positron polarisation measurement technique, have been pursued 
at our Institute^, ^, |^. Results are promising^, 0, and compete well with limits 
on physics beyond the Standard Model obtained from experiments at much higher 
energies. 

Typically in such (3 decay experiments, one is interested in measuring the longitu- 
dinal polarisation of positrons emitted parallel or antiparallel to the direction in which 
the decaying nucleus is polarised. The emitted positron is stopped in some specific 
medium where positronium is formed. The medium being placed in a strong magnetic 
fieldQ, the decay spectrum of positronium formed provides information on the po- 
larisation of the incoming positron. The analysis of actual experimental positronium 
spectra usually assumes that, at the location where positronium is formed, the ap- 
plied magnetic field is exactly parallel to the direction in which the positron is emitted. 
Since, strictly speaking, such an assumption can never be correct in practice, the more 
general situation has to be considered, namely when the direction of the magnetic field 
is arbitrary with respect to the positron momentum. However, the present author has 
been unable to find in the literature |^ ^, ||, |l^ a detailed discussion of this point, 
hence this note addressing the problem specifically. 

In sect. 2, a brief review of positronium hyperfine structure in the absence of any 
external electromagnetic field is presented, with also the purpose of specifying our no- 
tation. Sect. 3 develops the discussion of the positronium ground state and its hyperfine 
structure in the presence of an arbitrary static magnetic field. Sects. 4 and 5 consider 
lifetimes and time evolution of hyperfine populations, respectively, while sect. 6 ends 
with some conclusions. 



^ Since positrons and electrons have opposite electric charges, an external electric field docs not 
afifict positronium states nor their intrinsic decay rates (however, "pick-up" processes of external 
electrons or other "quenching" effects can affect these ratesg). Nevertheless, an external electric field 
can indeed improve the positronium formation rate in a given medium. 
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2 Hyperfine Structure of Positronium 



As is well known, the IS positronium ground state is in fact split into two hyperfine 
levels of total spin S = and 5" = 1. Their energy difference is due to spin-spin 
interactions between the positron and the electron, to relativistic corrections to their 
kinetic energies and to virtual pair annihilation in the S = 1 channel^. 

A non relativistic representation of the associated wave functions is sufficient for our 
purposes, as well as being justified. Accordingly, wave functions separate into space 
and spin components, with the space component being identical for both hyperfine 
states. Namely, the S = singlet or parapositronium state is given by a wave function 
of the form0 

I 0,0 >= iPir) -L [ 1+ > I- > - I- > 1+ > ] . (1) 

Here, ip{r) is the space wave function of the 15 ground state, while the second factor in 
the r.h.s. of this expression is the spin wave function. The convention used throughout 
is that in the product of two spin ket vectors, the first always stands for the electron 
spin component whereas the second stands for the positron spin component. These 
components are taken with respect, say to the z axis of some reference frame. Later 
on, this axis will of course correspond to the direction in which the f3 decay positron is 
emitted. It is also assumed that the space wave function ilj{r) is properly normalised, 

POO 

Att / drr^l ^{r) 1^= 1 , (2) 
Jo 

and that the basis vectors |+ > and |— > in spin space both for the electron and for 
the positron are normalised in the usual way, namely 

< +1+ > = 1 = < -|- > , < +|- > = = < -|+ > . (3) 

Consequently, the | 0, > state in ([l|) is also of norm 1. 

Similarly, the three components m = ±1,0 of the S = 1 triplet orthopositronium 
state are given by the normalised state vectors 

I 1, m = 1 >= ip{r) !+>!+>, 

I 1, m = -1 >= 'ijj{r) |- > |- > , (4) 

I 1, m = >= ■?/'(r) ^ [ 1+ > |- > + |- > 1+ > ] . 
v2 

The hyperfine states in (0) and @ are eigenstates of the total Hamiltonian Hq of the 
system in the absence of an external magnetic field. In the non relativistic limit, this 
Hamiltonian is comprised of the ordinary Schrodinger type Hamiltonian -ff(c) — which, 
apart from the usual kinetic term, only includes the Coulomb interaction between the 

^The latter effect is absent in the 5* = channel due to charge conjugation. A photon has C ~ — 1, 
whereas the S = and 5 = 1 states have C = +1 and C = — 1, respectively. 

^Note that a compactified notation for state vectors is used throughout, whereby their isotropic 
radial dependence is not displayed explicitly. 
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two oppositely charged particles — to which diverse relativistic corrections are added. 
The latter contributions correspond to spin-spin interactions between the electron and 
the positron, to relativistic corrections to their kinetic energies, to virtual pair anni- 
hilation effects for orthopositronium states and to further QED radiative corrections. 
Restricting for a moment the discussion to the non relativistic purely Coulomb Hamil- 
tonian, the singlet and triplet states above are degenerate eigenstates of -ff(c) with 
energyQ 

(5) 



(C) 



-"-2 2 

--a mc 
4 



-6.803 eV . 



In fact, in this purely Coulomb limit, the space wave function is simply 



jca 



-1/2 



-r/a 



with the positronium radius 



2h 



amc 



1.0584 A . 



(6) 



(7) 



Even though the complete wave function ilj{r) in (|l|) departs from the simple radial 
dependence in (^), the parameter a in (|^) gives a measure of the spatial extension of 
the positronium ground state. 

Relativistic effects just mentioned lift the singlet-triplet degeneracy^] according to 
the eigenvalues 



i/o I 0, >= ^0 I 0, > , Ho \ l,m = ±1,0 >= El \ l,m = ±1, > 
with energies [B| 



En 



-7 +" 
4 



1 5 
4 ~ 64 



El 



-7 +" 
4 



1 _ 5 1 
12 ~ 64 ^ 4 



1 1 
a mc 



2 2 

a mc 



(8) 

(9) 

(10) 



and the hyperfine difference |jT3|, 



AE = El- Eo 



|-f (f + 21n2) + In a-^ + 0{a^) 



14 2 

■^a mc 



lAl X 10"^ eV 



(11) 



Except for their first term corresponding to the purely Coulomb contribution (|^), 
the different contributions of order in (^ and ([l0|) are as follows. Both in Eq and 
in El, the first such contribution is that of the spin-spin interaction energy, while the 

'^Throughout, a, m and c of course stand for the fine structure constant, the electron and positron 
mass and the speed of fight, respectively. Numerical values for these parameters are from Ref. 

^In the purely Coulomb situation of (^, the singlet- triplet degeneracy is due|l^] to a dynami- 
cal SO{4:) symmetry explicitly broken by effects now considered. Nevertheless, the S — 1 triplet 
states m = ±1, remain degenerate since, in the absence of external electromagnetic fields, the total 
Hamiltonian Hq is invariant under rotations in space. 
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second — common to both expressions — is that due to relativistic corrections to the 
total kinetic energy. Finally, the third term of order in Ei follows from virtual pair 
annihilation in the orthopositronium triplet state. 

To conclude, let us consider positronium lifetimes. Due to charge conjugation prop- 
erties, the singlet state can only decay into an even number of photons and the triplet 
states into an odd number (beginning of course with three photons). Therefore in very 
good approximation, parapositronium decays predominantly into two photons and or- 
thopositronium into three photons, since compared to each of these processes, the rate 
for any further photon pair emission is suppressed each time by an additional power of 
a^. Hence, only 27 and 87 decay processes are considered in this note, and these two 
channels are assumed to encompass all possible decay modes of positronium. 

The associated lifetimes, including radiative corrections, have been computed within 
QED[|I5|, |T6|. For the singlet state, one has the 27 decay rate 

Xs =|«'^ [l-^{5-^) + laHna-' + 0{a' 

= (0.125209 X 10-9 s)-i = 7.98665 x 10^ s^^ . 
Similarly for triplet states, their 87 decay rate is 



(12) 



9tt 



a^^i-K^-d) 1 - (10.266 ±0.008)f-i-2 



(18) 



142.074 X 10-9 s)-i = 7.08859 x 10*^ s" 



Note the rather large ratio 



A5 
At 



1184.695 . 



(14) 



3 Coupling to a Magnetic Field 

Let us now consider positronium states formed in the presence of some external static 
magnetic field B = {B^, By, B^). For all practical purposes, certainly always reahsed 
in actual experimental conditions, it will be assumed that whenever B{x) might have 
a non vanishing gradient, this gradient is nevertheless negligible on the scale of the 
spatial extension of the positronium state, namely 

a ^ 1 , for all i = x,y, z . (15) 

I B I 

Here, a is the positronium radius of (^. Effectively, one may then consider the magnetic 
field B not only to be static but also to be constant, which is thus the situation to be 
assumed in the analysis developed in this note. Furthermore, the field B is not taken to 
be necessarily parallel to the z axis with respect to which spin eigenstates were defined 
in the previous section, since in practical applications, the latter axis is often defined 
by the positron momentum instead. 



4 



The presence of the magnetic field B induces an additional interaction term in the 
total Hamiltonian for the positronium system. The total Hamiltonian H now includes 
the previous Hamiltonian Hq with its lowest energy eigenstates in (0) and @, to which 
the magnetic coupling to magnetic moments is added, namely 



• B — jl^ ■ B . 



(16) 



Here, and /2_|_ are the electron and positron magnetic moments, respectively. In 
terms of their spin operators (T_/2 and (?+/2, respectively, we haveQ 



The magnetic dipole moment /i is given by 



1^ = 9 



eh 



with the gyromagnetic factor |18 



2m 

MM 



a 

2h + — + 

1 27r 



|C(3)-3C(2)ln2 + -C(2) + ^ 



- ) +0 



Therefore, the magnetic energy contribution to the total Hamiltonian reads 



(17) 



;i8) 



(19) 



(20) 



To complete this list of notations, it turns out that the parameter setting the 
physical scale of magnetic fields in the present system, is the combination 



2fi 



1 



AE 3.628575 Tesla 
leading to the definition of the positive quantity 



(21) 



X 



2fi 
'AE 



B 



B 



3.63 Tesla 



(22) 



In addition, it proves convenient[] to introduce the following combinations of the B^ 
and By components of the magnetic field B, 



B, 



Bx + iBy 
~V2~ 



5_ 



Bx — iBy 
~V2~ 



(23) 



Given the total Hamiltonian 



H = Ho + H 



IB 



(24) 



^Here, a are the usual Pauli matrices defining the spin 1/2 representation of the (double covering 
SU{2)) of the rotation group 50(3) in three dimensions. 

^The actual reason why these definitions are convenient is the fact that it is the spin 1/2 represen- 
tation of the three dimensional rotation group which appears throughout. 
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it is now a simple matter to proceed diagonalising it for its lowest energy states with 
spherical symmetry, namely in the sector of positronium 15" hyperfine states of the 
previous section. First, one easily finds 

I 0, >= I 0, > +fxB, I 1, > I 1, 1 > +fxB+ | 1, -1 > , (25) 

if I 1,0 >= El I 1,0 > I 0,0 > , (26) 
I 1, 1 >= El I 1, 1 > -fxB+ I 0, > , (27) 
H\l,-1 >= El I 1, -1 > +At5_ I 0, > . (28) 

Given these results, eigenstates of H and their eigenvalues can be derived after 
some work. The corresponding four eigenstates are denoted | ips' >, \ ipT' > and 
I ip± >. In the limit of vanishing magnetic field B, these states reduce — possibly up 
to some phase — to the singlet | 0,0 >, the triplet | 1,0 > and the triplet | 1,±1 > 
states, respectively, hence the notation. In particular, the | ipS' > and | ipT' > states 
are referred to as the "pseudo-singlet" and "pseudo-triplet" states, respectively. 

The pseudo-singlet state is given by 



I ^ps, >= -ij (1 + x^)-^/\ViT^ + 1)-'/' X {{VTT^ + 1) I 0, > 

-HE, I 1, > I 1, 1 > -115+ I 1, -1 >} , 

with the eigenvalue 



(29) 



Eci = — AE 
* 2 



Vl + x2 - ll +Eo . (30) 



The pseudo-triplet state is 

I >= ^ (1 + x^y^^WTT^ - 1)-'/' X {{VTT^ - 1) I 0, > + 

+^B, I 1, > I 1, 1 > I 1, -1 >} , 

with the eigenvalue 

Et' = +-AE 
2 

Finally, the remaining two states are 



(31) 



Vl + x^ + 1 +Eo. (32) 



, J_ B,±\B\ B- 111^ I J_ --B^l^l B+ I 1 _1 ^ 



(33) 



with the degenerate eigenvalue 

E± = AE + Eo = El . (34) 



Here, upper (resp. lower) signs in the r.h.s. of (|33|) correspond to the state | ip^ > 
(resp. \ ip- >). 
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By construction, the four states just given not only diagonalise the total Hamilto- 
nian H, but are also orthonormalised. Namely, these states are orthogonal by pairs and 
are normalised to 1. Of course, orthogonality is automatic for non degenerate states 
but not for degenerate ones; by construction, the states | ■?/'+ > and | > given above 
do indeed have a vanishing inner product. 

Before commenting on these expressions, it is useful to consider the limit in which 
the magnetic field B is parallel to the z axis, namely when 3^ = and By = 0. 
Of course, in this limit, the eigenvalues of the four states above remain unchanged, 
since their energies only depend on the variable x, i.e. the norm of the magnetic field. 
However, the states specified above, which diagonalise the total Hamiltonian H, then 
reduce to 

I >= COS0 I 0,0 > -sin0 I 1,0 > , (35) 

I >=r] { sine I 0,0 > +cos^ I 1,0 > } , (36) 

I V'i >= ie^^"^ \l,±V> ■ (37) 

In these expressions, rj = sign{Bz) is the sign of Bz, uj is some arbitrary phase whose 
value is dependent on the manner in which the limit B^ = 0, By = is takenQ, and 
the mixing angle 6 is defined by 

cos9 = ^ Jl + , ^ = , sin^ = r7^./l , ^ (38) 

Of course, these expressions coincide with those usually found in the literature]^, |T^, ^, 
in which case it is customary to take | il}± >=| 1, ±1 > and the z axis along the magnetic 
field S, namely ?7 = +1. In particular, note that in the limit of an infinite magnetic 
field -B, the pseudo-singlet and pseudo-triplet states in (^) and (^) further reduce to 
(in the notation of (|l])) 

I V'S' >= ^'^ V^(^) I >l +^ > , I V'T' >= '^{f) I +7] >| —f] > . (39) 



In order to comment on the physical significance of these results, let us first consider 
the case where the magnetic field is parallel to the z axis. Even though the vector B 
then explicitly breaks rotational invariance of the positronium system in vacuum, there 
still remains the symmetry of arbitrary rotations around the z axis. Consequently, the 
spin projection m on that axis still defines a good quantum number for positronium 
states (of vanishing angular momentum.). Therefore, the states | 1, 1 > and | 1,-1 > 
must remain eigenstates of the total Hamiltonian, whereas the other two states with 
m = 0, namely | 0,0 > and | 1,0 >, are now allowed to mix. Moreover, for the 
former two states with m = ±1, since the electron and positron have their spins then 
aligned and since their magnetic moments are equal in norm but opposite, the magnetic 
coupling energy Hb vanishes identically, leading for these two states to the same energy 
eigenvalue Ei as in the absence of any field B. Finally, for the remaining two states 

^Indeed, the coefficients of the states | 1, 1 > and | 1, — f > in ( ^3| ) are non analytic functions of 
and i?„. 
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with m = diagonalising the total Hamiltonian in the hmit^ where the magnetic 
couphng Hb is much larger than all other contributions to H , namely for magnetic 
fields whose magnitude is much larger than {AE/(2^) = 3.63 Tesla), the state of 
lowest (resp. highest) energy, i.e. \ ips' > (resp. | ipT' >), is the one for which the 
electron and positron magnetic moments are both aligned parallel (resp. antiparallel) 
to the magnetic field B, or equivalently the one for which the electron spin is aligned 
antiparallel (resp. parallel) to B and the positron spin is parallel (resp. antiparallel) 
to B. 

These properties — indeed, solely expected on physical grounds independently of any 
explicit calculation — are beautifully confirmed by the results above in the case where 
B^ = and By = 0. At this stage, it thus appears that an explicit calculation serves 
the purpose only of determining the mixing angle 6 in (RS) between the two m = 



states, and of deriving the energy eigenvalues Es' and Et' in (|30|) and (|3^), as functions 
of the magnetic field B = (0, 0, B^). 

Consider now the general situation when the direction of the magnetic field B is 
arbitrary with respect to the z axis. In fact, the results just discussed can be used 
in order to understand the general well. Indeed, the only difference between 

the two configurations is that the axis with respect to which the spin part of state 
vectors is expanded, is different. Hence, by an appropriate change of basis in the spin 
sector, effected through a rotation in the spin 1/2 representation, the eigenstates of 
the total Hamiltonian H in the arbitrary case {B^^By) ^ (0,0) can in principle be 
constructed from the expressions of these states when {B^^By) = (0,0). Therefore, 
since under such a rotation in spin space representations of spin and of spin 1 are 
invariant, only the states | l,m = 0, ±1 > can mix among themselves. Consequently, 
in the general case, both the pseudo-singlet and pseudo-triplet states should be given 
as some superposition of all four states | 0, > and | 1, m = 0, ±1 >, with in particular 
the coefficient of the | 0, > component independent of the components of the magnetic 
field but only dependent on its norm \ B \, whereas the remaining two states | ijj± > 
can only involve the three states | 1, m = 0, ±1 >. In addition, for all four eigenstates, 
the coefficients of the states | 1, m = 0, ±1 > must depend on all three components B^, 
By and B^ of the magnetic field. 

Indeed, these are features of the results in (PD|), (^TJ) and (|S3|), which therefore find 
their origin in the fact that spin and spin 1 representations are invariant under space 
rotations. However, only an explicit calculation — either along the lines just sketched or 
by direct diagonalisation of the Hamiltonian as done in this note — can determine the 
specific mixing coefficients defining each of the eigenstates of the total Hamiltonian H. 
Incidentally, note that the same argument of invariance under space rotations explains 
why eigenvalues of H remain unchanged when the magnetic field B is no longer parallel 
to the z axis, i.e. why these eigenvalues only depend on the norm | 5 | of the magnetic 
field. 



This situation was pointed out to the author by J. Dcutsch. 
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4 Positronium Lifetimes 



As a first application of results so far, let us compute decay rates for all four eigenstates 
of the total Hamiltonian H in the presence of a magnetic field B. Actually, such a 
calculation is rather straightforward, given the decay rates Xs and At of the singlet 
and triplet states | 0, > and | 1, m = 0, ±1 >, respectively. 

First, consider decay rates into two photons. Due to charge conjugation, only the 
I 0, > state has 27 decays. Therefore, the 27 decay rate of the pseudo-singlet state 
I ^5' > is 



x(27) 
A5, 



1 + 



Ac 



Ac cos 6 . 



For the pseudo-triplet state | ipT' >, we have 

1 



x(27) 

Arpi 



A, 



Ac sin^ e . 



(40) 



(41) 



Finally, the 27 decay rate of the remaining two states | il)± > is 

a1'^^ = . (42) 
Similarly, for the 87 decay rates of these states in the same order, one finds 

1 



x(37) 



X(37) 

Arpi 



1 

2 

1 
2 



At sin 9 



1 + 



At = At cos^ 9 



X, 



(43) 

(44) 
(45) 



In these expressions, the angle 9 is defined in (PB|). 

Finally, total decay rates — ignoring the much suppressed rates into four or more 
photons — are simply 



Ac' 



At' 



-{Xs + At) + - I A Xs - At 



2^^^ + ^^)-27rT^ 



= A5 cos^ 6* At sin^ 6^ 
{Xs — At) = Xs sin^ 6* + At cos^ 9 



A+ — Xn 



(46) 



(47) 
(48) 



As ought to be expected, these expressions only depend on the norm of the magnetic 
field 5, but not on its direction. Indeed, since the choice of axis with respect to which 
spin states are expanded does not affect the calculation of decay rates[^ that axis can 

^"indeed, positronium states can "remember" the direction and polarisation of the incoming positron 
only through the populations of the four Hamiltonian eigenstates (this is the topic of the next section). 
Decay rates are intrinsic properties of each of these states, and as such, are thus independent of any 
variable possibly affecting positronium formation. 
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always be taken along the magnetic field, in which case decay rates can only depend 
on I Bz I, namely the norm of the magnetic field. Incidentally, note that total 27 
and 37 decay rates — Xs and At, respectively — are independent of B — a consequence 
of unitarity. 

It is also easy to check that the statistical decay rates into two and three photons, 
i.e. the average of each of these rates over the four Hamiltonian eigenstates, are 
independent of the magnetic field. These statistical rates thus coincide with their values 
when B = 0, namely A5/4 and 3At/4 for two and three photon decays, respectively. 

Finally, let us remark that the difference between the total decay rates for the 
pseudo-triplet and | 4'± > states. 



At' — At — 2'''^'^ ~ '^'^^ 



VTT 



(As - At) sin^ 9 



(49) 



is a quantity always positive for all values of the magnetic field. 



5 Positronium Populations 

Let us now address the specific topic of this note; the formation of positronium in 
a medium placed in a magnetic field. In view of applications, the incoming positron 
is assumed to have a polarisation P+ (—1 < P+ < +1). This polarisation P+ is the 
expectation value of the positron spin projected onto its momentum, the latter vector 
thus also defining the z axis for spin quantisation from now on. Therefore, up to a 
physically irrelevant overall phase, the spin component of the incoming positron wave 
function is given by 

^^1 + ^+ I + > + I - > • (50) 

Here, 0+ is an arbitrary phase difference — thus possibly leading to observable physical 
effects — between the two spin components defining a positron state of polarisation P+. 

Similarly, it will be assumedQ that the positron capturing electron has a polarisation 
P_ (— 1 < P_ < +1) along the same z axis. Though in most practical applications, the 
positronium formation medium is at temperatures such that electrons are effectively 
not polarised, some experiments at very low temperatures are being planned, for which 
an investigation of possible effects due to electron polarisation in the applied magnetic 
field might therefore be found useful. Consequently, again up to a physically irrelevant 
overall phase, the spin part of the electron wave function is also of the form 

-^^TTT^ I + > + e^^-^Vl-P- I - > , (51) 

where 0_ is another arbitrary phase shift. 

Hence, at the moment of positronium formation, it is assumed that the spin com- 
ponent of the positronium state vector | = > is simply given by the tensor 

^^The interest of this possibility was pointed out to the author by F. Gimeno-Nogues. 
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product^ of the electron and positron spin vectors in (^Up and (0), while the space 
part of the state vector is of course the wave function ip{r) in ([^). In order to obtain 
the time evolution of the associated state vector, and thus also the time dependence of 
its decay products, the resulting wave function | = > of formed positronium has 
to be expanded in the basis of eigenstates of the total Hamiltonian H in the presence 
of the magnetic field B. This change of basis thus defines coefficients C5/, Ct' and C± 
such that 

\^,t = 0>= J2 Ca \^a> . (52) 
a=S',T',+,- 

Explicit expressions for these coefficients are given in Appendix A. Time evolution of 
the positronium state formed is then given by 

\i^,t>= Ca \'lpa> exp (-^Eat - ^Xat) , (53) 

with the quantities Ea and (a = S',T',+,-) defined in (jS^), (H, (H and (||), 
m and (H). 

Given these expressions, time evolution of each of the populations associated to 
each of the four states \ ipa > {cl = S', T', +, — ) is simply obtained as 

\Ca\' e-^"*, a = S',T',+,-. (54) 

Expressions for all observables of interest are then easily written down. For example, 
27 and 87 production rates are0, respectively, 

Ri2,)it)= E X^^'UCal' e-'^\ (55) 

a=S',T',+ - 

R(3.){t)= E A^)|ape-^»S (56) 

a=S',T',+ - 

while the total photon production rate — simply the sum of the latter two rates — is 
itself 

Rit) = >^a\Ca\' e-'^' . (57) 

a=S',T',+,- 

The rather lengthy expressions for the populations at t = 0, namely the coefficients 
I Ca P (a = S' ,T' , +, —), are given in Appendix B. Results probably more relevant 
at this point are the same coefficients averaged^ over the phase shifts 0_ and 
Indeed, one ought to expect that for most sources where the /3 decay process responsible 



^^Note that it is always possible to "rotate away" one of the two phases 0_ or (f)^ — but not both — 
by an appropriate rotation around the z axis. However, dependence on the cancelled phase then 
reappears through the and By components of the rotated magnetic field. 

-'^^Note that when actual experimental data are considered, the 2j production rate in (|55| ) should 
also include a "fast" component due to direct pair annihilation of incoming positrons with electrons 
of the positronium forming medium. 

-'^^This average does not amount to setting e*"^* = in the original expressions for the coefficients 
Cq, and can only be applied once the complete expressions of Appendix B for the coefficients | Ca P 
have been obtained. 
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for positron production is taking place, positron states with all possible phase shifts 
0+ are statistically populated, thereby justifying an average of the final positronium 
populations over 0+. Similarly, in the positronium forming medium, one should also 
expect that all electron states with different phase shifts 0_ are statistically populated, 
again justifying an average over the phase 0_. Under these assumptions^^, the averaged 
populations are given by 



C. 



S'T' 



[1 - P.P^ 



1 



2/i 



B. [P- - P4 



1 

+ 2 



It 



B 



-P^P. 



+ ; 



and 



In the r.h.s. 
(resp. 



I B I 



l±^P, 
\ B\ 



(58) 
(59) 



of these two expressions, the upper (resp. lower) sign refers to | Cs' |^ 
^ (resp. I C- P), respectively. Note that the sum of the four 

1, as ought be the case since, by 



Ct' n and I C+ |2 (resp. | C. 
phase averaged populations does indeed reduce to 
definition, the state \ ip,t = > is normalised to 1. 

In turn, phase averaged photon production rates R{2-y) 



(t), -R(3^)(t) and R(t) are 



K37) 

defined as in (^5]), (|56D and (0), of course involving now the phase averaged coefficients 
I Ca P (a = S",T', +,— ) just given. Clearly, these averaged photon rates are more 
readily amenable to experimental measurement than are the non phase averaged rates 
considered previously. 



6 Conclusions 

This note reports on the calculation of hyperfine positronium populations formed in the 
presence of an arbitrary magnetic field — whose gradient is assumed to be vanishingly 
small over the spatial extension of the positronium bound state — for arbitrary positron 
and electron polarisations. The analysis generalises results available!^, |10|, |[ in the 
literature in two respects. On the one hand, the magnetic field is not assumed to be 
necessarily aligned along the momentum of the incoming positron. On the other hand, 
allowing for possible electron polarisation effects enables the present results to be of 
application to positron polarimeters operated at very low temperatures. 

Expressions derived in this note provide the basic information required in any ex- 
perimental analysis of results obtained using positron polarimeters based on either 
time or energy distributions of positronium decay photons. As a simple but explicit 
illustration of relevance to current /? decay experiments]^, ^, 0], let us consider for 

^^Note that taking such averages is even more justifiable in the instance — often realised in practice — 
that the magnetic fields present in an experimental set-up possess an axial symmetry along the axis 
of incoming positrons. Indeed, as was noticed previously, either of the phases (f)- or (j)^ can always be 
"rotated away" by an appropriate rotation around the z axis, then also a symmetry transformation 
of the magnetic fields. 
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example the phase averaged photon production time spectrum R{t) when electrons in 
the positronium formation medium are not polarised (P_ = 0). Expressions derived in 
the previous section then lead to 



R{t) 



+ -rXi 



2fi 



2/i 



(60) 



However, since the pseudo-singlet decay rate A5/ is much larger than the pseudo-triplet 
one At', the pseudo-singlet contribution in the r.h.s. of ( |BUD becomes effectively neg- 
ligible after a nanosecond or so, leaving only the first two terms. Hence in effect, the 
photon time spectrum R{t) 



R{t > 1 ns) 
with the analysing power 



1 



Xtc 



-Xrt 



1, 

+ -At^ 



[1 - eP^ 



-\rp! t 



2fi 



(61) 



(62) 



provides the means of measuring positron polarisations. Note that, when compared 
to the usual result obtained for a magnetic field B assumed to be parallel to the 
incoming positron momentum, the sole effect of non vanishing components B^ and By 
in this simple example is to decrease the effective analysing power e multiplying the 
positron polarisation P+. Nevertheless, there certainly exist other instances where the 
effects of non vanishing components B^ and By have to be properly accounted for when 
analysing actual experimental data. The results of this note provide the basis for such 
investigations. 
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Appendix A 

This Appendix gives the coefficients of the hnear combination of the total Hamiltonian 
eigenstates defining the state vector of positronium formed in the presence of a magnetic 
field B (see (|5|)). 

For the pseudo-singlet and pseudo-triplet states, one has 



Cs'T' 



V2 



(i+x^)-^/^ (v/rT^±i)-i/2 



X 



2V2 



1 + p_ J I - p+ 



1 - P a/1 + e'' 



2V2 



1 + P_ ^Jl-P+ e'^+ + sJl-P- ^l + P+ e 



i(j> 



(Vl + a;2± 1) ^ 
2/i 



5. ± 



5+ T 



2 V V + 



5_ 



(63) 



where the upper (resp. lower) sign refers to the coefficient C5/ (resp. Ct')- 
Similarly, the coefficients of the remaining two states | ip± > are 



a 



± 



1 + P Jl - p. e''>'+ + Jl - PJl + P+ e 



i4>~ 



Bl + Bl 



B 



+ 



+ 



2V2 



1 + P_ a/1 + P 



P.± I B 



B, 



B 



Bl + Bl 



(64) 



+ 



2^2 



1 - P_ ./l - P, e^('^-+'^+) ^ 



-p,± 


B 


B^ 




B 


^Bl + 



Here again, the upper (resp. lower) sign refers to the coefficient C+ (resp. C_). 
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Appendix B 



This Appendix gives the populations of positronium states formed in the presence of 
a magnetic field, for electrons and positrons of initial polarisation P_ and P+, respec- 
tively. With the same conventions as to upper and lower signs as in Appendix A, the 
results are as follows. 

The I V'S'.T' > populations are 

X I 1(1 + x^Y/\^/lT^ ± 1) [1 - ^^-^+1 T 

T^(^rT^±i)|^s,[p_-p+]+ (%\\b_b^)p_p^t 



AE 



\AEJ 



1(V1 + X2±1)^1-P2 ^l_p2 _ ^^-) _ 



(65) 



4 \AEj 
2V2 \AE 

^ 2V2 \aeJ 



-icf)- 



I -Pi 



l_p2 



B_e''^+ + B+e-'^+ 
B_e''>'- + B+e^'^-' 



2fi 



VTT^±l)T^B,P 
2/i 



± 



Similarly, the | ip± > populations are 

\c^\'= l4^[i+p_p+]± l^[p_+pj+ 

' ' 41 S |2 ^ +J 41 S I ^ +^ 2 I S |2 



± 



4a/2 



1±4^P_ 



± -1-Jl-P^ 
4^2^ 



1 ± 



P I 

A 
I p 



B , 
P_ 
I P 



P I 
I P 



-40+ 



-»</>_ 



± 



+ 



B 



B 



P 



P 



(66) 



Note that in each of these two equations, the last three terms could be combined 
further into the product of two terms. However, results in the form given here are 
more amenable to the phase average discussed in the main text. Incidentally, it is 
straightforward to check that the sum of all four populations does indeed reduce to 1, 
since, by construction, the state \ ip,t = > is normalised to 1. 
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